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Pressure  versus  time  data  obtained  from  igniting  and  burning  propellants  in 
highly  loaded  closed  chambers  (loading  densities  between  0.1  and  0.3  g/cm3)  are 
being  used  to  evaluate  propellants  in  empirical  quality  and  production  control 
procedures.  Interpretations  obtained  from  such  production  control  devices  are 
hampered  by  uncertainties  concerning  ignition  delays,  dynamic  burning,  heat  loss-, 
propellant  surface  area  versus  distance  burned,  and  real  gas  effects.  In  this 
study,  attention  is  focused  on  devising  a closed  chamber  experiment  for  measuring 
dynamic  burning  rates  and  on  developing  a comprehensive  analytical  model  of  the 
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TMta  In  which  propellant  granules  are  ignited  and  Durned  in  closed 
chambers  are  use d empirically  as  part  of  routine  quality  and  production 
control  procedures  (to  determine  propellant  acceptability,  to  establish 
charge  weights,  arid  to  sake  laboratory  determinations  of  propellant  char- 
acteristics) .1-3  Closed  chamber  experiments  of  interest  in  this  study  use 
high  loading  of  nitrocellulose  based  propellant  (e.g.,  0.1  to  0.3  q/m*) 
which  produce  high  pressurization  rates  (dp/dt)MX  of  1"»*  MPa/sec  and 
greater)  and  high  pressure!  (e.g.,  100  to  .;00  MPa).  In  this  study,  we 
focused  our  attention  on  closed  chamber  experiments  a-,  a means  of  learning 
more  about  dynamic  bur nine  and  dynamic  chamlier  responses  that  accompany 
ignition  and  rapid  pressurization. 

The  reader  should  make  the  distinction  between  two  "<  ry  ditfmr.t 
situations  in  which  dynamic  Lurning  and  dynamic  chamber  lusimses  are 
of  interest: 

1.  Production  Control  - In  operational  systems  in  which  the  prenur- 
isation  processes  are  uniform  and  performance  is  r.ot  abnormal  1'/ 
sensitive  to  the  propellant  variations  that  oci.ur  dunna  produc*-icr, 
relative  changes  in  dp/dt  vs  p (uf rally  dunra  the  middle  40 

to  80%  of  the  pressure  range)  may  provule  sufficient  infcrnwtie ' 
to  establish  the  charge  weiqht  for  ea^h  production  propellant  1<  t . 
However,  further  improvements  of  muzzle  v#l<r it”  prediction:"  ha.  '*1 
on  closed  chamber  results  require  a m-.a  complete  understand1 1/" 
and  interpretation  oi  ;losed  chamber  ,rrccsf<. 

2.  Research  and  Development  - In  develop*- nr »'  ms  (e.g.,  where 

maximum  volumetric  loadinq  or  optimized  :inuli  j^-ition  inr.iterc 
are  sought),  careful  attention  must  be  i,ver.  to  the  dynamic 
chamber  and  burninq  responses  associates  with  ignition  and  initial 
pressurization,  tn  these  cases,  it  i'  ro  glean  igr.it  i.i-i, 

dynamic  burning,  and  chamber  dynamics  infcrm-iti'  r frem  sp-iciall", 
devised  closed  chamber  exiieriments.  Thus,  emphasis  is  given  to 
the  initial  portion  of  the  p vs  t trace. 

Recently,  attention  is  bcinq;  focused  or.  *.peciti<-  n tuat icrs  which  would 
benefit  from  improved  understanding  of  how  propellant  variations  affet 
closed  chamber  processes,  e.q.,  Ref.  1 is  directed  it  methods  of  predict m-j 
charge  weights  for  production  systems  and  Ref.  4 rr  j ited  nonuniform  prvrs  . e 
gradients  to  overpressure  failures  and  abnormal  muc'le  velocity  and  ranor 
dispersions. 

While  the  component  pre  censes  that  occur  in  cloned  chamhc-i  expel  im<  .t: 
are  reasonably  well  understood,  i.one  of  the  r-xictinq  analytical  models  trovide 
a comprehensive  explanation  of  the  interplays  amonq  these  component  j rorer,.-:es . 
For  example,  uncertainties  exist  concemino  the  contributions  of  dynamic  burn- 
ing, real  gas  effects  and  the  use  of  simple  equation*:  of  state,  heat  lenses 
from  the  gases  to  the  chamber,  variations  of  flame  end  chamber  temperatures, and 
nonuniform  ignition  and  surface  roqression  of  granulir  propellant  chorees. 


For  the  aoat  part,  previous  efforts  considered  the  middle  of  the  burning 
interval  to  avoid  complications  associated  with  the  above  stated  uncertain* 
ties.  Approaches  which  do  not  treat  adequately  the  ignition  and  early 
pressurization  phases  ignore  siuch  of  the  information  available  from  closed 
chamber  experiment*. 

By  understanding  the  dynasiic  responses  in  the  ignition  and  early  pressur- 
ization phases,  important  insights  into  the  instabilities  and  oscillations 
that  have  been  associated  with  Ml functions  in  operational  systems  (4)  can  be 
obtained.  In  particular,  the  goals  of  the  study  jure: 

1.  To  develop  ignition  and  dynamic  burning  rate  relationships  that  can 
re  tried  readily  in  interior  ballistic  codes. 

T.  To  jrov.de  an  analytical  basis  for  developing  an  improved  method  of 
measuring  dynamic  responses  of  propel  Lants  at  hiah  pressures  and 
high  pm.-unzation  rates,  and 

*.  To  develop  a more  complete  model  of  the  interactions  between  ignition, 
dynamic  burning,  chamber  pressurization,  heat  loss.  etc. 

DESCRIPTION  OF  CLOSED  CHAMBER 

closed  chamber:  are  usually  designed  to  burn  granular  propellant  charges. 
Several  characteristics  of  these  devices  complicate  the  data  interpretation : 

1.  Granular  propellants  ire  usually  ignited  by  an  unevenly  dispersed 
secondary  solid  proj'ellant  charge  and,  thus,  uncertainties  axist 
concerning  tin*  uniformity  of  ignition  and  flame  spreading. 

C . Heat  Ions  iron  the  combustion  gases  to  tne  chamber  walls. 

. Burning  suitor-  arti  versus  distance  burned  is  altered  by  non uni form 
ignition  p.r.1  tnus  r.o.t uni  form  surface  regression)  and  variations 
m tht  muJt  .-perforated  grain  geometry. 

4.  Ignition  ci»arg«.-.  usually  produce  gas  compositions  and  temperatures 
that  dufe’  Hiomficantly  from  those  of  the  main  charge. 

For  experimental  determinations  of  batch- to-batch  acceptability  during  propel- 
lant production,  the  above  four  compli cat ions  may  not  be  objectionable  since 
unacceptable  prop  »llur.t  can  l>e  isolated  and  charge  weight  determined  in  terms 
of  relative  changed  in  the  pressure  versus  timt  trace.  However,  the  complica- 
tions completely  mash  changes  in  burning  rate,  and  in  particular,  the  dynamic 
responses  that  occur  during  ignition  and  the  early  phases  of  chamber  pressur- 
ization. 

We  propose  a special  purpose  closed-chamber  and  propellant  configur- 
ation which  is  suited  tor  establishing  the  contributions  of  dynamic  burning 
and  chamber  responses.  As  shown  in  Fig.  1,  the  propellant  charge  consists 
of  thin  sheets  of  propellant  bonded  to  the  inner  wall  of  the  cylindrical 
chamber  and  to  tp,,.  outer  diameter  of  the  cylindrical  spacer.  For  this  simple 
geometry,  the  surface  area  vorvus  distance  burned  is  known  with  a relatively 
high  degree  of  certainty.  For  the  datum  case  configuration  of  Table  I,  the 
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surface  area  decreases  6.5%  as  the  propellant  burns.  Heat  loss  to  the 
chaaber  walls  does  not  appreciably  affec*;  chamber  pressure  and  temperature 
*^nc*  propellant  insulates  the  walls.  We  propose  to  greatly  reduce  the 
complication  associated  with  nonuniform  ignition  and  with  the  uncertainties 
of  the  gas  composition  by  using  an  igniter  charge  that  consists  of  ultra-thin 
shavings  (i.e.,  5 to  15ji  thick)  of  the  mai.<  propellant.  These  shavings  will 
loosely  fill  the  chamber  free  volume  and  will  be  ignited  by  an  electrically 
heated  nichrome  wire;  flame  will  propagate  quickly  through  the  bed  of  shavings 
and  provide  an  intense  heat  flux  to  the  propellant  surface.  The  rate  at 
which  the  propellant  shavings  burn  can  be  increased  by  pre-pressurizing  the 
chamber.  Heat  flux  from  the  ignition  charge  will  be  measured  using  an  outer 
wall  section  instrumented  with  thin  film  calorimeters.  The  chamber  will  be 
fitted  with  a blow-out  disk  to  limit  the  maximum  pressure  whenever  limited- 
***»9e  pressure  transducers  are  used  to  obtaii  very  accurate  pressure  measure- 
ments at  the  lower  pressures. 

Figure  lb  shows  an  alternate  configuration  in  which  the  propellant  con- 
figuration is  a hollow  cylinder  which  burns  on  the  inner  and  outer  cylindri- 
cal walls  and  the  ends.  The  advantage  of  this  configuration  is  that  the 
propellant  charge  is  easier  to  fabricate;  the  disadvantage  is  that  the 
chaaber  walls  are  fully  exposed  and,  thus,  the  heat  loss  will  be  greater. 

ANALYSIS 


The  chamber  pressure  and  mass  generation  responses  (indicated  in  Fig.lc) 
will  be  calculated  from  a mathematical  model  whose  primary  components  are 
coupled  through  numerical  solutions  to  the  energy  equation  for  the  propellant, 
the  Zeldovich  nonsteady  heat  feedback  equations,  the  energy  and  continuity 
equations  for  the  chamber,  and  a generalized  form  of  the  equation  of  state. 

Burning  Rate  Dynamics 

During  ignition  and  rapid  pressurization,  two  opposing  effects  in  the 
combustion  zone  cause  the  instantaneous  burning  rate  to  differ  greatly  from 
the  steady  state  burning  rate  corresponding  to  the  instantaneous  pressure, 
i.e.,  (1)  the  developing  temperature  profile  in  the  condensed  phase  of  the 

propellant  and  (2)  the  out-of'phase  blowing  effect  of  the  reacting  gases 
leaving  the  burning  surface.  The  contributions  of  these  effects  are  mojit 
easily  explained  in  terms  of  an  example  which  considers  a rapid  pressure 
increase  from  a low  pressure  condition.  At  the  lower  burning  rate  corres- 
ponding to  the  lower  pressure  prior  to  pressurization,  the  thermal  wave  is 
relatively  thick.  (Similarly,  slow  ignition  produces  a thermal  wave  thicker 
than  the  thermal  wave  corresponding  to  steady  state  burning.)  As  the  pres- 
sure increases,  the  propellant  surface  layer  is  in  effect  preheated  with 
respect  to  the  thermal  profiles  at  the  higher  pressures;  the  burning  rate  is 
enhanced  while  the  preheated  surface  layer  is  being  burned  away.  This  en- 
hanced burning  rate  causes  the  flame  zone  gases  to  have  an  enhanced  velocity 
and,  thus,  the  "blown"  flame  is  thicker  than  the  steady  state  flame  corres- 
ponding to  the  instantaneous  pressure.  This  blowing  effect  decreases  the 
temperature  gradient  at  the  surface  and  the  heat  feedback  from  the  flame. 

Thus,  the  burning  rate  of  preheated  propellant  is  moderated  by  the  blowing 
effect  when  rne  burning  rate  exceeds  the  equilibrium  burning  rate.  The 
interactions  resulting  from  the  rapid  changes  in  pressure  will  also  affect  the 
flame  temperature  and  chamber  gas  temperature.  A primary  objective  of  this 
study  is  to  elucidate  these  dynamic  processes. 
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In  this  development,  the  dynamic  burning  rate  calculations  employ  the 
Zeldovich  formulism  (5)  which  offers  the  important  advantage  of  considering 
the  burning  race  transients  without  having  to  understand  the  details  of  the 
flame  and  surface  zone  reactions.  Since  the  methodology  is  covered  thorough- 
ly in  Ref.  5,  only  the  application  is  discussed  in  this  paper. 

The  analysis  and  assumptions  are  discussed  as  they  apply  to  highly 
loaded  closed  chambers.  Because  of  the  wide  range  of  events  and  special 
purpose  propellants  that  may  be  considered,  the  assumptions  used  in  the 
analysis  should  be  reviewed  for  each  situation  considered.  The  standard 
assxsnptions  for  application  of  the  Zeldovich  nonsteady  heat  feedback  function 


1.  The  rate  processes  in  the  flame  zone  adjacent  to  the  surface  (i.e., 
the  fizz  zone)  and  m the  surface  reaction  zones  can  be  considered 
quasi-steady  in  the  sense  that  their  characteristic  times  are  short 
compared  to  the  transition  time  of  the  pressure  change.  Mo  such 
limitation  is  necessary  for  the  thermal  wave  in  the  condensed  phase. 

2.  No  kinetic  heat  release  occurs  in  the  condensed  phase  below  the 
surface  reaction  zone.  Although  the  surface  reactions  occur  in  a 
zone  of  finite  tmekress,  the  zone  is  sufficiently  thin  that  it  can 
be  considered  as  quasi-steady. 

3.  The  condensed  phas«  of  the  propellant  is  accurately  represented 
as  being  homogeneous  and  isotropic. 

4.  Propellant  combustion  zones  are  not  influenced  by  position  and 
external  forces  such  as  shear  forces  from  the  flowing  chamber  gases. 

The  justifications  for  assumption  1 closely  follow  the  arguments  pre- 
sented in  Ref.  5.  Thus  following  Ref.  5,  the  characteristic  times  for  the 
condensed  phase,  surface  zone,  and  gas  phase  aro: 

t « a /r2  (1) 

c c 

x = KT  /E  t - 0.1  t (2) 

s s s c c 


Tf  = [VcPf/(V;£°c)lTc  = °'2  X 10~  pTc 


For  the  datum  case  conditions: 


at  p = 100  atm 


= 500 


--  C.0003  sec  and  t-  = 0.02t  sec 


= 0. 00004 


- 0.1t 


Thus,  over  the  pressure  range  of  significant  dynamic  burning  effects,  the 
flame  zone  responds  very  quickly  when  compared  to  the  response  of  the  con- 
densed phase  temperature  profile.  Furthermore,  the  flame  zone  reacts  quickly 
relative  to  the  overall  pressurization  time  (typically,  0.001  to  0.008  sec). 
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It  ,f  , 


The  energy  equation  for  the  condensed  phase  (-•  < x < 0)  has  the  fol- 


lowing eigenvalue  dependence  on  Tc(0,t.)  * Tg  (t) 


pc  (3T  /at  + r (T  )3T  /3x]  * X 32T  /3x2 
c c c sc  c c 


where  under  conditions  of  full  ignition  the  burning  rate  r is  related  to 
the  surface  tesqperature  through  the  pyrolysis  law. 


r = Ag  exp(-Eg/RTs) 


Note  that  surface  temperature  is  dominated  by  the  heat  feedback  from  the 
flame.  The  initial  condition  for  Eq.  (4)  is. 


T (x)  * T, 


t = 0 


The  first  boundary  condition  is 


T = T„ 


x 


The  second  boundary  condition  is  a series  of  conditions: 


(1)  heatup  to  ignition 


For  T T 
s 


ze'  *3xJ, 


(Until  onset  of  surface  regression,  the  temperature  profile  is  calcu- 
lated by  a closed  form  heat  conduction  solution.) 


(2)  combined  heating  from  ignition  stimulus  and  from  gas  phase  and 
surface  reactions 


For  T > T 

s — ze 


3t1 

and  t < ti  , X— I ■=  q (t)  + Xc$  [r,p] 
3 J x=0  3 


(3)  burning  dur  .ng  conditions  of  full  ignition  after  the  igniter  burns 
out 


For  T > T and  t > t.  , a—  = > 4>[r,p] 
s — ze  — ig'  3x1  - cY 

3 ' x^O 


where  $ is  the  Zeldovich  heat  feedback  function  developed  in 
Ref.  5. 


The  result  of  the  foregoing  analysis  of  the  propellant  ignition  and 
burning  is  one  .nonlinear,  second  order  partial  differential  equation, 

[Eq.  (4)]  for  T{x  t) , with  a series  of  boundary  conditions. 


An  expression  for  the  dynamic  flame  temperature  may  be  obtained  under 
the  following  conditions:  (1)  the  product  composition  in  the  flame  zone  and 


H 
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the  heats  of  reaction  in  the  surface  reaction  zone  and  in  the  flame  zone 
are  not  affected  by  the  dynamic  conditions,  and  (2)  the  specific  heats 
in  the  surface  reaction  zone  and  the  flame  are  equal.  It  follows  from 
t * (r/ac) (Ts  - Tq)  that  the  energy  balance  for  both  steady  and  unsteady 
states  is: 

rpc£hf  - V + *0ccc<Tf  - V <U) 

Therefore  by  considering  a reference  condition,  the  dynam:c  flame  temperature 


£ °=cp,E  rJ 


—It  . - t 1 + t. 

j cp,fLs'r,:£  SJ  1,1 


Chamber  Pressurization  Dynamic? 


Hie  analysis  and  assumptions  are  discussed  as  they  apply  to  closed 
chambers  in  which  the  combustion  gases  are  well  stirred,  i.e.,  thermal 
gradients  within  the  chamber  are  not  considered  to  be  important.  The  as- 
sumptions for  the  closed  chambers  shown  in  Fig.  1 are: 

1.  Temperature  and  pressure  variations  can  be  accounted  for  by  a 
lumped  parameter  treatment. 

2.  The  flame  zone  gases  are  fully  reacted  when  they  enter  the  control 
volisne  which  is  adjacent  to  burning  propellant  surface. 

3.  Propellant  surface  ignites  uniformly. 

The  condition  of  negligible  total  temperature  and  pressure  gradients 
of  assumption  1 is  more  justified  for  our  small  free  volume  chamber  (defined 
by  Fig.  1 and  Table  I)  with  a distributed  propellant  charge  than  for  a 
large  free  volume  chamber  with  the  propellant  charge  at  one  end.  In  the 
former  case,  the  combustion  gases  leaving  the  burning  propellant  surface 
are  distributed  along  the  length  of  the  chamber  and  have  velocity  components 
(normal  to  the  axis)  which  promote  rapid  mixing. 

Assanption  2 is  valid  for  nitrocellulose  propellants  at  sufficiently 
high  pressure  that  the  dark  zone  thickness  is  small  compared  to  the  dimen- 
sions of  the  free  volizne  surrounding  the  burning  propellant.  Establishing 
the  validity  of  this  assumption  will  have  to  await  analysis  of  experimental 
data,  since  during  the  early  phases  of  pressurization  of  high  volumetric 
loading  chambers  reactions  in  the  chamber  may  be  significant. 

The  validity  of  assumption  3 depends  on  devising  the  experiment  so 
that  the  interval  of  flame  spreading  through  the  propellant  shavings  is 
short  compared  to  the  total  burning  time  of  the  main  propellant  charge. 

The  derivation  of  the  chamber  equations  was  carried  out  without  speci- 
fying the  gas  law.  During  the  numerical  solution,  gas  properties  are  cal- 
culated as  part  of  the  finite  difference  solution  (either  from  tabular 
data  or  from  an  equation  of  state  and  thermodynamic  relationships)  for  each 
set  of  calculated  pressure,  p and  temperature,  T^.  The  initial 
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conditions  are 


f 


pch  ~ pinit 


T = T 
ch  init 


(13) 


The  mass  continuity  equation  for  the  free  volume  in  the  chamber  (see  chamber 
control  volume  in  Fig.  la)  is 

d<(,chV/dt  * "b  + "ig  1141 

”h«”  % ’ PcW1' 

The  chamber  volume  increase  during  burning  is 

dV  . /dt  * rA.  + m.  /p 
ch  b ig  rc 

Following  ignition,  the  burning  surface  area  is  a prescribed  function  of 
propellant  geometry  in  terms  of  distance  burned,  xfa. 

The  energy  equation  in  the  chamber  is 

dt  = pv  [mb <hf  “ e)  + mig(hig  “ e)  Closs] 
ch 


(16) 


(16) 


For  M-6  propellant  a simpler  form  of  the  energy  equation  can  be  used  with 
very  little  loss  in  accuracy  since  hf,  cv  and  dh/3T)p  vary  slightly 
(see  Table  II)  and,  as  shown  in  Fig.  2,  the  Noble-Abel  equation  of  state 


is  an  excellent  approximation.  Accordingly,  de  - c dT,  and 


dT  m 1 

dt  c pV 

v ch 


(m.  (c  T-  - c T)  + m.  (c  T.  - c T)  - Q 
b p f v ig  p ig  v 


loss 


(17) 


ttie  generalized  equation  of  state  p = p/Tw(p,T)  is  expressed  in  differential 
form  (and  then  combined  algebraically  with  Eq.  14): 


dp  1 dp  p dT  p du 

dt  uT  dt  _2  dt  _ 2 dt 
uT  Tu 


(18) 


In  the  numerical  solution,  tu(p,T)  may  be  input  in  tabular  form  and  particu- 
lar values  of  u(p,T)  are  found  by  a two  point  interpretation. 


The  heat  loss  from  the  chamber  gases  to  the  chamber  walls  are  approximated 
by  treating  the  chamber  wall  as  several  thermally  thick  surfaces  that  are 
heated  by  a combination  of  radiation  and  convection 


Q 


loss 


+ CA 


1-2 


(IS) 
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The  wall  temperatures  are  calculated  by  integral  solutions  to  the  heat  con- 
duction equation  [for  discussion  of  Eq.  (20)  see  Ref.  6] 


T 

w,n 


4a 


3A 


VTch  - V 


(T  - T-)  (2T  . - T - T ) 
w 0 ch  w 0 


dt 


(20) 


Equation  (20)  was  obtained  by  substituting  a third  order  temperature  profile 
into  the  heat  conduction  equation  with  a convective  boundary  condition  and 
setting  up  the  indicated  integration. 

As  the  igniter  gases  pass  over  the  unignited  propellant,  they  give  up 
heat.  As  a result,  the  igniter  gas  tes^erature  as  it  mixes  with  the  chamber 
gases  is  considered  to  be 


T.  = T.  , - A.  q.  /(c  m.  ) (21) 

ig  ig,ref  d ig  p ig 

The  primary  result  of  the  foregoing  analysis  of  the  chamber  and  propel- 
lant burning  conditions  is  a coupled  set  of  three  nonlinear,  first  order 
ordinary  differential  equations  (Eq.  j.4,  15,  and  17)  for  pCn#  Tch*  and 
Vch  and  one  nonlinear,  second  order 'partial  differential  equation,  i.e., 
the  energy  equation  for  the  condensed  phase,  Eq.  (4).  The  energy  equation 
for  the  condensed  phase  along  with  its  boundary  conditions  was  solved  using 
the  methods  of  explicit  finite  differences.  Equations  (14),  (15),  and  (17) 
were  solved  by  a predictor-corrector  technique  employing  a variable  time 
step. 


CALCULATED  RESULTS 


A series  of  parametric  studies  were  carried  out  for  the  closed  chamber 
configuration  shown  in  Fig.  la  loaded  with  M-6  (a  single  base  nitrocellulose 
propellant).  The  datum  case  combustor  conditions  listed  in  Table  I,  were 
arrived  at  after  using  the  analytical  model  to  investigate  a range  of  con- 
figurations. The  M-6  condensed  phase  properties  selected  for  the  study  are 
representative  of  typical  nitrocellulose  based  propellant  properties;  a 
consistent  set  of  property  measurements  for  M-6  has  yet  to  be  accomplished. 
The  combustion  gas  properties  were  calculated  using  the  techniques  of  Ref. 

7.  As  indicated  in  Table  II,  the  ranges  of  M*,  3h/3T)p,  cv  and  Tv  are 
surprisingly  small  and,  thus,  mean  property  simplifications  to  the  equations 
are  soundly  based.  A similar  series  of  thermochemical  calculations  for  a 
double  base  propellant  revealed  that  the  ranges  of  these  same  variables  are 
approximately  twice  as  great.  Figure  2 indicates  that  the  large  departures 
from  ideal  gas  theory  and  that  in  the  range  of  interest, single  real  gas 
laws  (e.g.,  the  Noble-Abel  equation)  can  adequately  describe  the  gas  proper- 
ties. 


Burning  rate,  chamber  temperature,  and  pressure  during  the  operating 
sequence  of  the  datun  case  are  shown  on  Fig.  3a.  After  the  propellant  ignites, 
the  heat  flux  from  the  ignition  charge  produces  a burning  rate  overshoot 
since  it  is  the  same  order  of  magnitude  as  heat  feedback  from  the  propellant 
flame.  The  second  overshoot  in  burning  rate  is  the  increase  in  burning  rate 
produced  by  the  rapid  increase  in  chamber  pressure.  Note  that  burning  rate 


and  chamber  temperature  dynamics  are  very  prominent  during  the  first  one- 
third  of  the  burning  interval.  The  short  tail-off  time  of  the  igniter  charge 
(i.e.,  O.S  msec)  produces  the  sharp  demarkation  between  the  burning  rate 
overshoot  produced  by  the  igniter  and  the  burning  rate  overshoot  during  the 
subsequent  rapid  pressurization n extending  the  tail-off  time  to  1.0  msec 
merges  the  two  overshoots  and  there  is  no  clearly  distinguishable  demarka- 
tion between  the  overdriven  burning  rate  during  ignition  and  ditring  rapid 
pressurization  (see  the  comparison  on  Fig.  3b). 


The  most  widely  used  result  from  closed  chamber  experiments  is  dp/dt  at 
prescribed  pressures.  Figure  4 shows  the  strong  influence  that  loading  den- 
sity has  on  dp/dt.  Figure  5 illustrates  the  application  of  a widely  used 
parameter#  relative  quickness  (RQ) , to  relate  the  dp/dt  vs  p results  of 
closed  diember  experiments.  One  of  the  more  widely  used  definitions  of 
relative  quickness  (RQ)  states  that  dp/dt  be  recorded  at  four  equally 
spaced  values  of  pressure  in  the  middle  portion  of  the  pressure-versus-time 
interval. (2)  The  dp/dt's  are  then  compared  to  datum  cube  or  standard  lot 
values  at  the  four  pressures  in  the  following  manner. 


RQ 


4 

I l 

4 n*l  dp/dt) 


dp/dt 


datum 


(22) 


For  the  four  pressures  indicated  in  Fig.  5b  (i.e..  Pi  M 500,  P2  * iC‘00, 

P3  * 1500,  and  P4  » 2000),  RQ  is  referred  to  as  RQo. 5-2.0* 

Figure  5 shows  how  changes  in  burning  rate  exponent  (while  keeping  the 
burning  rate  at  68  atm  constant)  affect  dp/dt  a.n.1  RQo. 5-2.0*  Wh«n  consid- 
ering the  effects  perturbing  n,  it  must  be  realized  that  the  values  of 
RQ  are  very  dependent  on  the  pressure  at  which  burning  rate  is  held  constant 
(i.e.,  the  pivot  point,  p,  on  a In  r versus  In  p plot).  Indeed,  as 
the  pivot  point  pressure  is  increased  toward  the  mean  value  of  p^,  p£»  P3 
and  P4,  3RQ/3r.  becomes  very  small. 

Figure  6 illustrates  that  the  Eq.  (22)  definition  of  RQo. 5-2.0 
totally  blind  to  variations  that  greatly  affect  ignition  and  the  initial 
pressurization  phases.  Even  though  decreasing  the  igniter  heat  flux  and 
mass  significantly  alters  the  burning  rate  and  chamber  temperature  responses 
before  p * 300  atm,  after  p = 500  atm  is  achieved  the  changes  in  dp/dt  vs  p 
are  not  appreciable. 

Figure  7 and  Table  III  summarize  the  extent  that  changing  propellant 
parameters,  the  closed  chamber  configuration,  and  the  model  affect  the 
ignition  and  pressurization  parameters.  Based  on  the  RQo. 5-2.0  definition, 
variations  in  the  burning  rate  coefficient  and  exponent  have  the  greatest 
influence  on  RQ.  To  illustrate  how  dynamic  response  data  can  be  obtained 
from  closed  chamber  experiments,  Table  III  includes  RQo. 05-0.2  calculations 
based  on  pressures  at  50,  100,  150,  and  200  atm.  Note  that  these  values  of 
RQo. 05-0. 2 strongly  reflect  changes  in  dynamic  response. 


Figure  8 is  a Zeldovich  map  which  shows  burning  rate  as  a function  of 
heat  feedback  to  the  condensed  phase  of  the  propellant.  Following  the  Zeldo- 
vich formalism  nonsteady  burning  rates  can  be  related  to  steady  state  burning 


rata  changes  that  correspond  to  chan ges  in  the  j ropellant's  initial  tempera- 
ture. A detailed  discussion  of  this  relationship  is  given  in  Ref.  5.  The 
locus  of  burning  rates  illustrate  the  relationships  between  burning  rate  over- 
shoot and  increases  in  initial  tesperature.  For  example,  the  r/r«q  * 2 over- 
shoot following  ignition  (see  Fig.  2)  corresponds  to  the  increase  in  steady 
stats  burning  rets  obtained  by  .increasing  the  initial  tenperature  of  the  pro- 
pellant fraa  298  to  approximately  480  K.  As  the  dynamic  effects  diminish,  the 
locus  of  blaming  rates  merges  with  the  burning  rates  corresponding  to  T0  ■ 298  K. 

* 1 

Several  recent  efforts  have  attempted  to  develop  empirical  dynamic  burning 
relationships  in  terms  of  dp/dt  end  other  instantaneous  parameters.  One  such 
relationship  originated  with  von  Elbe  (6)  and  has  been  extended  by  others  (a.g.. 

Rtfs.  1 and  9).  To  data  none  of  these  empirical  relationships  successfully  com- 
pered with  the  more  complete  dynamic  burning  rate  models. (5,10)  Two  of  the  sever- 
al reasons  for  this  are  (1)  if  the  method  ct  lur.  h is  used,  the.  essential  fea- 
tures cf  the  condensed  phase  temperature  history  are  lost  since  th«-  computational 
scheme  uses  a series  of  steady  state  temg.eratu*e  profiler  in  the  -endensed  phase, 
end  (2)  the  method  does  not  consider  the  dynamic  coupling  between  the  heat  feed- 
back from  the  flame  rone  and  the  proj  client  surface  conditions  (i.e.,  the  rela- 
tionship indicated  in  Eq.  10) . Two  examples  of  the  lack  of  effectiveness  of  the 
empirical  models  are  shown  in  Fig.  9.  If  the  empirical  relationships  are  to  bm 
useful,  an  approximate  one-to-one  relationship  between  r/r^  and  some  dp/dt 
parameter  must  be  found.  The  left-hand  figure  shows  that  input  parameters  that 
influence  the  dynamic  response  greatly  alter  the  relationship  between  r/r*q  end 
(<>c/r|q)  (d  In  p/dt) . Proponents  of  the  empirical  dynamic  relationship  may 
counter  that  since  changing  op  changes  the  propellant,  the  dynamic  burning  < 

relationship  will  change  also.  However,  in  the  right-hand  figure,  the  propellent 
parameters  are  the  same  for  all  three  cases,  but  iqnition  and  pressurization 
responses  have  been  altered.  Note  that  r/raci  is  very  dependent  on  the  partic- 
ular ignition  and  pressurization  response  and  is  not  specified  (even  approximate- 
ly) by  the  parameter  (ac/r£q) In  p/dt) . Attempts  to  correlate  other  dp/dt 
parameters  were  unsuccessful.  Thus  the  burden  of  demonstrating  that  <:  i>articu- 
ler  empirical  dynamic  burnir.q  rate  equation  is  a good  approximate  is  on  the 
usern  of  the  equation.  In  the  meantime,  there  is  -o  iustr f jc«t ion  ter  usina 
empirical  dynamic  burning  rite  equa*  ions  to  s-t*  if-  i . burring  rate 

responses- . 

CONCLUSIONS 

An  analytical  model  ha-  been  developed  that  considers  the  complete  thermo-  f 

dynamic  sequence  in  the  closed  chamber  configuration  that  ir  being  proceed  for  i 

future  studies  of  dynamic  n-.t-onses.  The  model  tro<r,s  fir  Le.ond  previous  efforts 
to  simulate  closed  chambers  smev  it  combines  for  tne  r'-.rst  time  the  following  . ! 

features: 

1.  Heat  up  to  ignition. 

2.  Dynamic  burr.ing  rate  and  flame  trrpc»'ature  Luseu  on  transient  devel- 
opment of  the  propellant  temperature  ircf-.le. 

3.  A generalized  treatment  of  real  gas  efforts. 

4.  Heat  loss  to  the  chamber  walls. 

5.  Chamber  wall  temperature  calculated  is  a function  of  time. 

6.  The  chamber  temperature  based  on  i solution  to  the  energy  equation. 
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tawd  on  the  analytical  results,  tha  operability  and  tha  Table  I dimen- 
aiona  for  tha  closed  chaabar  shown  in  Pig.  la  appear  to  ha  raaliatic.  Ef- 
forts are  proceeding  to  obtain  data  using  the  device  and  to  explain  the  data 
in  terns  of  the  model. 

Present  BMthods  of  Measuring  relative  quickness  exclude  the  contribution 
of  ignition  and  early  stage  pressurization  events  which  are  known  to  be  the 
precursors  to  serious  overpressures  and  velocity  dispersions  produced  by 
improperly  designed  large  caliber  propelling  chargee.  Accordingly,  whenever 
dynamic  response  information  is  sought,  careful  attention  must  be  given  to 
the  early  portion  of  the  p versus  t trace.  The  routine  quality  control 
type  of  HQ's  are  totally  blind  to  dynamic  response  irregularities  since 
they  use  data  from  the  well  behaved  middle  portion  of  the  dp/dt  versus  p 
trace. 


For  systems  with  smooth,  well  behaved  pressure  versus  time  characteris- 
tics, the  routine  quality  control  type  of  RQ  (based  on  the  middle  portion 
of  the  dp/dt  versus  p trace)  is  of  value  in  establishing  charge  weight. 
Under  such  conditions,  this  study  shows  that  variations  in  the  coefficient 
and  exponent  of  the  burning  rate  equation  influence  RQ  most  strongly  and 
that  the  assumptions  of  no  heat  loss,  Noble-Abel  equation  of  state,  no  dyn- 
amic burnino.  and  ignoring  the  details  of  the  ignition  processes  have  little 
influence  on  the  calculated  RQ  values  taken  from  the  middle  of  the  dp/dt 
versus  p trace. 
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Nomenclature 


burning  rate  law  coefficient  for  empirical  relation  apn 

2 


propellant  burning  surface 

pre-exponential  factor  in  pyrolysis  law 

overall  interchange  factor  for  radiation 
between  gas  and  region  being  heated 

specific  heat 

specific  heat  of  combustion  gases  at 
constant  pressure 

specific  heat  of  combustion  gases  at 
constant  volume 

activation  energy  in  pyrolysis  law 
enthalpy 

convective  heat  transfer  coefficient 

net  heat  release  during  propellant  combustion 

mass  flow  rate 

average  molecular  weight 

exponent  of  apn  burning  rate  law 

pressure 

heat  flux  from  ignition  stimulus 
universal  gas  constant 

relative  quickness  (see  discussion  of  Eq.  22) 

burning  rate 

time 

temperature 

initial  temperature  of  propellant 
threshold  surface  temperature  for  first  flame 
chamber  volume 
distance 

thermal  diffusivity 
thermal  conductivity 
density 

temperature  gradient  at  interface  between 
surface  reaction  zone  and  non-reacting 
condensed  phase 

Stef an- Boltzmann  constant 

temperature  sensitivity  of  burning  rate  at 

constant  pressure  (3ln  r/3TQ) 

covolume  in  Noble  Abel  equation  of  state 


cm 
cm/sec 


cal/g-X 

cal/g-X 

cal/g-X 

cal/g-mole 

cal/g 

2 

cal/cm  -sec-X 
cal/g 
g/sec 
g/g-mole 

atm>  MPa 

cal/cm2-sec 

cal/g-mole-X 

cm/sec 

sec 

K 

K 

K 

3 

cm  f 
cm 

cm2/sec 

cal/cm-K-sec 

g/cm3 


X/  cm 


K-1 

cm3/g 


-13- 


fl 

i 

I 

-4 


* *• 


'ts&axesz&auk  -- 


'^'Ti 


T 

characteristic  time 

see 

Tb 

thickness  of  propellant  burned 

cm 

y 

main  and  igniter  propellant  loading  density  in 
chamber 

g/'o 

Subscripts 

b 

burning  propellant 

c 

condensed  phase 

ch 

chamber  conditions 

«q 

steady  state  condition  which  corresponds  to 
instantaneous  pressure 

f 

flame  condition,  flame  zone 

ideal 

perfect  gas 

init 

initial;  condition  of  a controlled  environment 

ig 

igniter 

P 

constant  pressure 

prop 

propellant 

ref 

reference  conditions:  p = 68.08  atm  (6.90  MPa), 

N 

O 

s 

surface  reaction  zone 

V 

constant  volume 

w 

wall 

ze 

temperature  of  flame  development 

0 

ambient  condition  of  unburned  propellant 
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298  K 


Table  X 

Combustor  Conditions  For  Datta  Cat* 


loading  Density,  p * 0.2  g/ca 3 
Grain  Configuration!  rod  and  tub* 

Length  * 3.00  as 

Chaaber  ID  ■ 3.000  eta 

Rod  00  ■ 1.411  ca 

Neb  Thickness  - 0.050  cm 

Chaaber  Volume:  14.47  to  16.45  cm3 

Burning  Area:  43.05  to  40.26  ca2 

Igniter  Burning  Tiae:  1.0  msec  plus  0.50  Mac  tail-off 

Igniter  Charge  Mass  * 0.0187  g 

Heat  Flux  froa  Igniter  * 200  cal/ca2-sec 

Propellant  Mass  > 3.284  g 

Exposed  Hall  Surface  >11.0  ca2 

Chaaber  Material:  Steel 


Table  II 


Properties  for  M-6  Propellant  as  Used  in  Datum  Case 
COMBUSTION  GASES:  Range:  300-3300  at*  and  T * 2400  to  2800K 


M^  « 23.28  g/g-aole 
3h/3T)p  - 0.444  cal/g-K 
c >0.349  cal/g-K 


±0.086  (±  0.4%) 

±0.012  (±  2.7%) \ RANGE 
±0.017  (±  4.9%) 


pv/nRT  in  tabular  form  ranges  froa  1.033  tc  1.402 

T > 2821  K ranges  ±5K(±  0.2%)  as  loading  ranges  froa  0.1  to 
V 0.3  g/cm3 


CONDENSED  PHASE: 

r " rref(p/pref)n  cn/sec 
where  n * 0.674 

rre£  * 1.448  cm/sec  at  68.08  atm  (6.90  MPa)  and  298K 

0 - 0.0035  K”1 

P 

Eg  * 15,000  cal/g-moln 

T > 623K  at  1.448  ca/sec 
s 

X « 0.00073  cal/crs-sec-K  (±  0.00007) 

p > 1.58  g/cm?  (±  0.01) 

c « 0.37  cal/g-K  (±  0.03) 

T * 500K 
ze 
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Table  III 


' 1 

! 


Sensitivity  of  RQ  to  Datum  Case  Parameters 


Variable  Changed 

Fractional 

Change 

Ax/x 

RC0. 5-2.0 

ARQ 

Ax/x 

RC!0.05-0 

n 

+0.074 

1.15 

2.04 

1.036 

n 

-0.074 

0.87 

1.78 

0.965 

n 

+0.020 

1.038 

1.924 

1.C10 

n 

-0.021 

0.963 

1.860 

0.991 

a 

+0.020 

1.019 

0.988 

1.020 

a 

-0.020 

0.979 

1.023 

0.978 

U#  Loading  density 

+0.050 

1.060 

1.187 

1.106 

U , Loading  density 

-0.050 

0.941 

1.173 

0.948 

0 

+0.143 

1.001 

0.005 

1.048 

P 

a 

-0.143 

0.999 

0.004 

0.989 

P 

n»ig.  Igniter  mass 

-0.25 

1.001 

-0.003 

1.005 

qig,  Igniter  heat  flux 

-0.25 

0.999 

0.000 

0.932 

No  heat  loss  to 
end  wall 

— 

1.015 

— 

1.019 

Noble-Abel  Eqn. 
rather  thant>>(p,t) 

— 

0.998 

— 

0.994 

Increase  igniter 
tail-off  0.005  sec 

0.995 

M 

1.366 

T0 

+0.02 

1.027 

1.351 

1.040 

c,  Specific  heat 

-0.05 

1.000 

-0.007 

1.038 

X , Thermal  conductivity 

+0.05 

1.000 

0.006 

0.98-* 

(a) 

(b) 

ROD-AND-TUBE  PROPELLANT  CON- 

HOLLOW TUBE 

PROFELLANT  CON- 

FIGURATION (SIDE  WALLS  NOT 

FIGURATION 

FIDE  WALLS  ARE 

EXPOSED)  ■” 

EXPOSED) 

TUBEWALL 

(c)  CHAMBER  AND  PROPELLANT  PROCESSES 

y>  y 

k smfat  m.'tiu-m'VWM'ik  'e. . 

i 1 . *■'•!:  / 


j (^VTig  j: 

j j— - i “b^t  ^ 

, qlOS8  1 _ ^ 

PROPELLANT  T 


■RF 


Toi|i 


CONDENSED  PHASE 
SURFACE  'ZONE 
FLAME  ZONE 


TEMPERATURE  GRADIENT 
IS  NONSTEADY 


ENLARGED  SECTION 
OF  PROPELLANT 


Fig.  1 Closed  chamber  configurations  along  with  chamb 

and  propellant  processes  considered  in  mathemativ-d 
formulation. 


COMPRESSIBILITY  FACTOR,  pv,nRT 


PRESSURE,  MPa 
100  200  300 


M-b  PROPELLANT 


2400  K- 
2600  -v 
2800 


CALCULATED  BY 
BRL'S  BLAKE 
THERMOCHEMISTRY  CODE 

(RtF.  7) 


1000 


2000 


3000 


4000 


PRESSURE,  ATM 

Departure  of  M-6  propellant  combustion 
gases  from  perfect  qas  theoty  in  pressure 
and  temperature  range  of  closed  chamber 
experiments . (Note:  the  above  plot  the 

Noble-Atel  equation  of  .;tate  produces  a 
straight  line  with  a sope  of  hMy^/RT,  where 
b is  the  covolume.  Thus,  since  on  the  alo. e 
plot  the  isotherms  car  te  approximated  by 
such  a straight  line,  the  Noble-Abel  equation 
describes  the  gases  in  the  range  of  interest.) 


0 


PRESSURE,  ATM 

Fig.  2b  Isochoric  flame  temperatures  and 

molecular  weights  for  M-6  propellant 


PERIOD  or  IMPORTANT 
NONSTEADY  BURNING 
INTERACTIONS 


PERIOD  OF  RQ 
MEASUREMENTS 
USED  IN  PRODUCTION 
CONTROL 


If 

v max, ideal 

| IGNITION 

/ p . . = 2039  ATM 

^max, ideal 

= 206.6  MPa 

TIME  FROM  ONSET  OF  IGNITER,  MSEC 

Fig.  3a  Ignition  and  pressurization  transients  of  datum 
case  showing  burning  rate  overshoots  resulting 
from  igniter  and  pressurization. 
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— IGNITER  BURW  hijt  IN  1.5  MSEC 

IGNITER  BURN  OUT  IN  2.5  MSEC 


IGNITER 
sy/  :i/vSS  FLUX 


TIME  FROM  ONSET  OF  IGNITER,  MSEC 

3b  Extending  the  action  time  of  the  igniter  beyond 
the  time  of  the  burning  rate  overshoot  associated 
with  ignition  preheating  causes  the  two  burning 
rate  overshoots  shown  on  Pig.  3a  to  merge. 


r 

f 

1 

i 

0 

PRESSURE,  MPa 
100 

200 

* 

1 

dp/dt  x 10  ATM/bii; 


"i«j.  5a  Increasing  relative  cjirknoss  with  increasing 

burning  rate  exponent.  (Four  pressures  used  in 
relative  quickness;  calculation  are  indicated.) 


Fig.  5b  Increasing  pressurization  rate  with  increasing 
burning  rate  exponent. 


DATUM 

IGNITER  MASS  DECREASED  25% 
IGNITER  HEAT  FLUX  DECREASED  2i% 


Lp,ref 


/// 


* 

DATUM 
m.  -4 


7/.. 


IGNITER 
/yMASS  FLUX 


'0.5-2 


DATUM 
mig  -25% 

*ig  -25% 

^^'*p/,pmax,  ideal 

= 0.999  for  q^  -25% 
= 1.001  mig  -25% 


TIME  FROM  ONSET  OF  IGNITER,  MSEC 

6 Variations  in  igniter  produce  marked  changes 
in  the  ignition  transients  but  the  shape  of 
the  latter  portion  of  the  p - t curve,  and 
thus  relative  quickness,  is  not  affected. 
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RELATIVE 


NORMALIZED  BURNING  RATE 


LOCUS  OF  BURNING  RATE  /I 
- * 68  atm  = 6.90  MPa  / / / 

,ref  = 298  K / V 

ref  ' 1‘448  cn/3cc  / ,// 

/[/  //  F// 

/ A/'srf  / p 

'/A's'js' 

/l'  S 272  ATM 
//.  y (27.6  MPa) 


68  atm  =6.90  MPa 


*////}/  1™  ATM  (13.8  MPa) 


68  ATM  (6.9  MPa) 


M-n  DATUM  CASE 


37  ATM  (.1.72  MPa) 


Xc(dT/dx)s,c/\.(dT/dx)Si.r(,f 
NORMALIZED  HEAT  FLUX  TO  CONDENSED  PHASE 


8 Zeldovich  map  showing  traverse  of  datum  case 


H 
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f INCREASE 
. IGNITER 
SHUT  OFF 
BY  0.5 
MSEC 

DATUM 


LOADING 
-0.15  G/CM - 


